Spin phonon induced colinear order and magnetization plateaus in triangular and 
kagome antiferromagnets. Applications to CuFe02. 
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Coupling between spin and lattice degrees of freedom are important in geometrically frustrated 
magnets where they can lead to degeneracy lifting and novel orders. We show that moderate spin- 
lattice couplings in triangular and Kagome antiferromagnets can induce complex colinear magnetic 
orders. When classical Heisenberg spins on the triangular lattice are coupled to Einstein phonons, 
a rich variety of phases emerge, including the experimentally observed four sublattice state and the 
five sublattice l/5th plateau state seen in the magneto-electric material CuFe02. In addition we 
predict magnetization plateaus at 1/3, 3/7, 1/2, 3/5 and 5/7 at these couplings. Strong spin-lattice 
couplings induce a striped colinear state, seen in Qf-NaFe02 and MnBr2. On the Kagome lattice, 
moderate spin-lattice couplings induce colinear order, but an extensive degeneracy remains. 



Frustrated magnets, in which classical ground states 
have multiple accidental degeneracies, have been at the 
focus or renewed attention. While a large body of theo- 
retical work has focused on the lifting of this degeneracy 
by thermal or quantum fluctuations ('order by disorder' 
[l|) in real materials, alternate mechanisms might domi- 
nate. One mechanism that is always present is the cou- 
pling of magnetism to the lattice (spin-phonon coupling) . 
Exploring the ramifications of such couplings is still in its 
infancy and leads to new theoretical problems with con- 
sequences for real materials. Finally, the subject is po- 
tentially important for future applications - spin phonon 
couplings could give rise to multiferroic systems, where 
magnetic order is coupled to ferroelectricity 

The role of lattice distortions in promoting valence 
bond physics (spin Pierels effect) is well known and con- 
sequences for spin 1/2 frustrated quantum magnets have 
been studied [3, 01- In the opposite large spin (semi- 
classical) limit, the role of lattice distortions in selecting 
the colinear ground states amongst the vast number of 
degenerate spin configurations of the pyrochlore magnet 
has been emphasized |5|- Further selection of a unique co- 
linear ground by various spin phonon interactions have 
also been proposed [1, 0] . 

Here we will study the effect of spin-phonon interac- 
tions on the classical ground states of the triangular and 
Kagome magnets. In contrast to earlier applications, 
all ground states of the antiferromagnets on these lat- 
tices are non-colinear (120°) configurations. Hence a fi- 
nite strength of spin phonon coupling is required to in- 
duce colinear states. We believe this to be the origin 
of the puzzling colinear magnetism seen in several trian- 
gular lattice materials, such as CuFe02[l], a-NaFe02[i], 
and MnBr9[loj. In these materials with S = 5/2 filled 
shell moments, magnetic anisotropy is likely to be (and 
in some cases known to be) very small, and cannot be 
invoked to explain the colinear order. Furthermore, we 
study the precise pattern of the colinear magnetic order 
induced within the Einstein Site Phonon (ESP) model 



0, which builds in an energy cost for displacing a mag- 
netic atom from its ideal lattice position, and is parame- 
terized by a single coupling constant. Analytic arguments 
combined with numerical calculations (simulated anneal- 
ing) are used to establish the zero temperature phase 
diagram of the classical spin models - which is briefiy 
summarized below. On the triangular lattice, increasing 
the spin-phonon coupling leads to the phase diagram in 
Figure 1, where the four sublattice zig-zag state (Estate) 
is the first colinear state selected. At large couplings a 
four sublattice parallel stripe state {S state) of up and 
down spins is obtained. The phase diagram in a magnetic 
field is remarkably complex, with magnetization plateaus 
at l/5th, l/3rd, 3/7th, 3/5th, 5/7th, and 1/2 occur in 
the parameter range where the Z state appears, and also 
l/9th and 7/15 at larger spin-phonon couplings. In par- 
ticular, the four sublattice Z state is typically converted 
into a l/5th plateau with a five sublattice spin structure 
on increasing the field. 

These are precisely the two colinear states observed 
in CuFe02, [8| an extensively studied triangular mag- 
net, where the complex four and five sublattice colinear 
structures described above was seen. This was previ- 
ously rather mysterious since the only available models 
which captured such orders were Ising Hamiltonians with 
large and very specific second and third neighbor inter- 
actions (J2 = O.45J1, J3 = 0.75Ji)[lH (or more compli- 
cated models, see e.g. Plumerfl^) which are unnatural 
since the magnetic susceptibility is nearly isotropic in the 
paramagnetic state. In contrast, our isotropic spin phon 
model involves a single parameter - the spin phonon cou- 
pling. Predictions for higher magnetization plateaus and 
other states for this material are made. The magnetically 
induced electrical polarization observed in this system is 
however not captured by our simple model, pointing to 
the role of other interactions which we briefiy discuss. 
The S state is precisely the spin arrangement seen at low 
temperatures in a— NaFe02 and MnBr2. Predictions for 
magnetic plateaus in these materials are also made. 
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On the Kagome lattice too we find that beyond a crit- 
ical coupling, colinear ground states are obtained, but in 
contrast to the triangular lattice, the manifold of these 
states has an extensive entropy. Only at larger couplings 
is there a transition into a unique ground state. This may 
be of relevance to the recently studied Kagome staircase 
compound, Mn3V2 08[l3'|. If the spin phonon mechanism 
is as widely relevant as some of our results suggest, the- 
oretical studies will have to move beyond their current 
focus on purely rigid lattice models. 

Spin Phonon Model: The minimal spin phonon cou- 
pling arises from the dependence of the exchange cou- 
pling on separation between the magnetic ions J(r). If 
sites i and j approach each other by Uij, then the spin 
phonon Hamiltonian then is: 



H = J 



E 

<y> 



(1 - aijUij)Si ■ Sj + iilatticc({ui}) (1) 



where a — J~^dJ/dr, -ffiatticc is the phonon Hamil- 
tonian and Ui is the displacement of site i. Hence, 
Uij = {ui — Uj) ■ Cij is the relative displacement pro- 
jected on the bond ij {cij is the unit vector from site 
i to j). Two types of phonon Hamiltonians have been 
proposed. First, the bond phonon model of Penc et al.[5| 
where iJiatticc = (-^/2) I]<ij> "fj- This model assumes 
that bond displacements are independent, which may be 
an oversimplification for many lattices. Integrating out 
the phonon degree of freedom here generates just the bi- 
quadratic term —bJJ2<ij> i^i ' ^j)'^- (The phonon dy- 
namics can be neglected if its frequency is much larger 
than magnetic energy scales). While for sufficiently large 
b, colinear states result, the model is still highly degener- 
ate (Ising antiferromagnet ground states on the triangu- 
lar or Kagome magnets). The selection by quantum fluc- 
tuations will be discussed elsewhere Instead, here 
we turn to the second model, the ESP model of Bergman 
et al.0 



iilatticc = {K/2) ^ , 



(2) 



which in contrast to the bond phonon model, respects 
the inevitable correlations between bond length arising 
from the underlying lattice structure. The ESP model 
describes a dispersionless optical phonon branch. More 
realistic phonon model informed by the detailed crys- 
tal structure of the material and acoustic phonons could 
generate longer-range effective spin interactions that are 
hard to deal with. In the interests of simplicity and gen- 
erality we will restrict attention to the ESP model. In- 
tegrating out the lattice displacement Ui results in the 
effective spin Hamiltonian: 



H. 



ESP 



= J 



J2 S.-S,-cJ2f^ 



<i]> 



(3) 



where c = J / {2K) is a positive dimensionless coupling 
and we have defined the dimensionless 'force' on site i as 
= Eje neighbors of ri^i ' Sj)eij, the sum in this defini- 
tion is over the (six) neighbors of site i. The spin-phonon 
interaction seeks to maximize the force which will 
result in gaining the maximum spin interaction energy. 
Note, the second term in Eqn. [3] generates interactions 
involving three adjacent spins e^j • ejk{Si ■ Sj){Sj ■ Sk)- 

Triangular Lattice: In the following we will con- 
sider a single triangular lattice sheet governed by the 
Hamiltonian in Eqn. ([3]) i.e. with nearest neighbor an- 
tiferromagnetic interactions and additional interactions 
generated by the spin-phonon term. We focus on the 
zero temperature phase diagram, as a function of the sin- 
gle parameter c and subsequently in an applied magnetic 
field. Since we are primarily interested in large spin (eg. 
S — 5/2), we focus on classical spins, where we can write 
Si = Sfii, where hi is a unit vector. This leads to an 
extended classical Heisenberg model of unit vectors with 
rescaled couphngs J = S^J and c = S^c . If we choose 
to restrict to Ising states. Si — <TiS, this effective spin 
Hamiltonian simplifies to an Ising model with nearest-, 
second- and third-neighbor coupling, J — cJ, cJ and c J 
respectively. However, in this model the second and third 
neighbor couplings are constrained to be strictly equal. 
This might be a rationalization for the large second- and 
third-neighbor couplings used in previous Ising models 

n|. 



Consider the ground state for classical spins on raising 
c. While at c = the the regular 120° pattern of 0(3) 
spins on the triangular lattice is realized, this is expected 
to survive to finite c as well. The ground state energy 
per site is Eq/J = —3/2, and since the force vanishes in 
this state, it is independent of c. While a full numerical 
solution is required (and provided below) for the phase 
diagram of this model, we can make some plausible ana- 
lytic arguments which will be confirmed by the numerics. 
Clearly, colinear states are preferred for large c since they 
give rise to the maximum force. However, near the phase 
boundary with the 120° states, the exchange J will pre- 
sumably be important, and hence we restrict attention 
to those colinear states that best satisfy J. These are 
nothing but the ground states of the triangular lattice 
Ising antiferromagnet (TLIA) with nearest neighbor ex- 
change, with two (one) up and one (two) down spins per 
triangle. We can then ask, which configuration within 
this manifold optimizes the force term? This question 
can be rigorously answered - it is the Z state. Using the 
dimer representation of the TLIA states, where a dimer 
is drawn orthogonal to each unsatisfied bond, and leads 
to a hard core dimer configuration on the honeycomb 
lattice, we see that the 'force' on site i is determined 
by the dimer configurations on the hexagon surround- 
ing site i. By checking the possible dimer configurations 
on a hexagon, we see that the force |_Fi| =2 if there is 
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FIG. 1: (Color online) zero field phase diagram of the tri- 
angular lattice ESP model. On increasing the spin phonon 
coupling c, the 120° state is followed first by the zig-zag Z 
state, then an 8 sublattice state and finally the stripe S state. 
Accidental degeneracies (thick red ticks) only occur at the 
phase boundaries c = 1/6, 1/2. Solid (hollow) blue (red) cir- 
cles are up (down) spins. 



one dimer in the hexagon or two not-opposite dimers; 
otherwise \Fi\ = 0. The ground state maximizes 
Hence it must have two dimers in every hexagon (since on 
average there are two dimers per hexagon and a configu- 
ration with a one dimer hexagon implies also a hexagon 
with three dimers, which experiences no force). Combin- 
ing this with the condition that the two dimers cannot 
be on opposite sides leads us to the unique zigzag state 
Z, with a four sublattice unit cell as shown in Fig. [1] and 
ground state energy per site Eq/J — —I — 4:C. 

The full phase diagram is obtained using simulated an- 
nealing on lattices with periodic boundary condition and 
various sizes up to 10 x 10, and choosing the state with 
lowest energy per site. Simulations on each size were 
done by an exponential annealing schedule from /3J = 0.1 
with a random initial state to (3 J = 1000, with a total 
of 20000^-^40000 sweeps, the whole process was repeated 
10 times to ensure stability of results. Site-update with 
Metropolis dynamics was used. While the algorithm does 
not guarantee convergence to the ground state we nev- 
ertheless believe an accurate picture emerges since all 
analytic expectations have been met, and we have not 
been able to guess ground states with better energies. 

The numerically obtained phase diagram in zero field is 
shown in Fig. [T] At c = 1/8 the 120° state is replaced by 
the four sublattice .Estate, with alternating zig-zags of up 
and down spins and ground state energy per spin Eq/J — 

— 1 — 4c. This is stable till c = 1/6, where an 8 sublattice 
state {Eq/J — —10c) takes over. Beyond c = 1/2, the 
S state of up-up-down-down stripes {Eq/J = 1 — 12c) is 
found and persists to large couplings. At the transition 
points c=l/6 and c=l/8, there are additional accidental 
degeneracies, and quantum effects could be important in 
resolving theseH^. 

The phase diagram in a magnetic field — 

— J^ihnl is remarkably rich. For small c where the 
120° state is realized, a highly degenerate set of states 
[l^ well known for the Heisenberg triangular antiferro- 
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FIG. 2: T = phase diagram of the triangular lattice ESP 
model. The vertical axis is the normalized field h — H/J. 
The fraction / label magnetization plateaus, while c— / are 
canted (non-colinear) states deriving from them. Also, cZ is 
the canted Z state while the 3 — sub states at small c also 
occur in the pure Heisenberg model in a field. 



magnet are obtained. Since they all have vanishing force 
contributions, the spin phonon interaction does not split 
this degeneracy. At larger values of c , the simulation 
shows a plethora of plateau states, which we briefly dis- 
cuss here and leave details to 15| . Interestingly the 1/5- 
plateau with the pattern observed in CuFe02 occur for 
a wide range of parameter c. For the parameter inter- 
val 0.14 < c < 1/6 our model shows both the zigzag 
Z ground state at zero field and the 1 / 5-plateau in mag- 
netic field, as in CuFe02. Other prominent plateaus that 
occur in the range of c where the Z state appears are the 
3/7th and 5/7th states with 7 site unit cell, a 1/2 magne- 
tization plateau with an 8 site unit cell and two distinct 
3/5th plateaus with 5 sites per unit cell. There is also a 
small region of 1/3 plateau, with a 12 site unit cell. The 
evolution of a plateaus state with increasing field can 
proceed in two ways. Either it can jump directly to an- 
other plateau, or undergo a canting transition, where the 
field direction is not parallel to the staggered moment. 
For example, the zero field Z state in the fully isotropic 
model undergoes a spin-flop transition immediately on 
applying a fleld, and the staggered moment is orthogo- 
nal to the field and the induced uniform magnetization. 
The evolution of the 1/5 plateau state on increasing the 
field however is continuous, with a gradual tilting of the 
staggered component away from the field. This phase 
boundary can be calculated analytically and agrees very 
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well with the simulations. Such canted states are labeled 
c — / in the figure (where / refers to the plateau they de- 
rive from) Such states are of course absent in Ising model 
studies ll| . Other plateaus occur for larger c, which will 



be discussed in detail elsewhere [15|. Here we note that 
the 1 /9th plateau that extends all the way down to zero 
field occurs because of the accidental degeneracy at the 
point c = 1/2 which includes finite magnetization config- 
urations, with a maximum magnetization of l/9th, which 
is selected by the field. Amusingly, the most obvious 1/3 
plateau expected for a triangular lattice, consisting of 
up, up and down spins on the three sublattices, does not 
occur (the 1/3 plateau at high fields involves a 12 spins). 

CuFe02 and other materials: In CuFe02 the 4 sub- 
lattice Z state is observed, which persists in a field upto 
B < 6 Tesla. At higher fields B > 14Tesla, the 5 sublat- 
tice l/5th magnetization plateau is observed. We note 
that both these states occur in our spin phonon model 
when 0.14 < c < 1/6. To estimate the spin phonon 
couphng c — c?J jlK in CuFe02, we use J = 39Kelvin 
(from the measured Weiss constant [l3]) and estimate 
a - 7 and if - 10,000Kelvin [J, which gives c - 0.1 
which is in the right ball park. While an isotropic spin 
model with magnetic order cannot have a magnetization 
plateau centered at zero field, even a very small magnetic 
anisotropy (e.g. an easy axis anisotropy — i'X^i ^Id ^^"^ 
produce the observed zero magnetization plateau, since 
the plateau width Ai3 is readily seen to be Ai? c>c V DJ. 
A 1% anisotropy Dj J produces the right plateau width 



17| . The magnetization profile as a function of field at 



c = 0.15 with a 1% easy axis anisotropy is shown in Fig. 
[3] (the field scale J/g^iBS is ~10Tesla per unit from the 
estimated value of J). The higher field magnetization 
plateaus and structures are predicted for future experi- 
ments on single crystals. Existing pulsed field measure- 
ments on powder samples reveal a sequence of anomalies 
at different fields up to full polarization, but the magneti- 
zation plateaus and structures remain to be conclusively 
identified [ll]. Lastly we note that the ferroelectric phase 
with incommensurate spiral order observed experimen- 
tally in CuFe02 in the field range 7 < B < 14Tesla[i3| 
is not produced here, indicating the importance of other 
lattice couplings eg. to the oxygen atoms mediating the 
superexchange interaction. The up up down down stripe 
pattern for 1/2 < c has been observed as the ground 
state for some materials with trian gula r lattice structure 
such as Q;-NaFe02,t2j and MnBr2, [lO|. Finding the pre- 
dicted l/9th and 1/3 plateaus in these materials would 
be a check of the spin phonon origin of these states. 

Kagome Lattice: The Einstein site phonon model on 
the Kagome is virtually identical to the triangular case, 
except that the lower symmetry in this case (lack of site 
centered 60° rotation symmetry) allows for anisotropic 
confining potential on the atoms. For simplicity, we as- 
sume an isotropic confining potential, but the main re- 
sults are independent of this assumption. 



1 

0.8 
0.6 
0.4 
0.2 



— I 1 1~ 

D, = 0.01,c = 0.15 



iiiiiiiiiiiii 




1 2 3 4 jj5 6 7 



10 



FIG. 3: (Color online) Predicted magnetization (m) curve for 
CuFe02 with c = 0.15 and 1% anisotropy. Plateaus occur at 
m = 0, 1/5, 1/3, 3/7, 1/2, 3/5, 5/7 and of course at m = 1. 
One unit of field is ~ 10 Tesla. Inset: the m = 1/5 state. 




FIG. 4: (Color online) zero field phase diagram of the ESP 
model on kagome lattice. Extensive degeneracy (marked red) 
persists into 1/12 < c < 1/6 where colinear states occur - 
representative configurations are shown. 



Simulated annealing was applied to this model with 
similar settings as the triangular case. The zero-field 
phase diagram is presented in Fig. [31 For small c 
we still get the ground states of the pure Heisenberg 
model, which are known to be extensively degenerate. 
For c > 1/12 we get colinear states, but in the range 
1/12 < c < 1/6 an extensive degeneracy remains [20j. 
Even more interestingly the zero-field ground states can 
have arbitrary magnetization ranging from —1/9 to 1/9 
per site. Therefore, in this zero temperature classical 
model, applying a small field will immediately induce a 
1/9-magnetization-plateau state. We expect that ther- 
mal and/or quantum fluctuation can lift this accidental 
degeneracy which is left for future work [l5| . Further in- 
creasing c beyond 1/6 pushes the system into a unique 
colinear states (see Fig. S]). 

We have shown that spin lattice couplings can induce a 
rich variety of complex colinear orders on the triangular 
and Kagome lattices. The moderate spin-phonon cou- 
plings strengths required make this a viable mechanism 
to explain similar ordering patterns seen in materials like 
CuFe02 . Our results underline the need to go beyond the 



current focus on purely rigid lattice spin models. Future 
work will study the effect of quantum fluctuations and 
valence bond states that naturally arise from disordering 
these colinear configurations. Including other atomic dis- 
placements might account for interesting magnetoelectric 
phenomena, and guide the search for multiferroic mate- 
rials. We acknowledge support from LBNL DOE-504108 
and useful discussions with K. Damle and R. Cava. 
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